Electric field control of valence tautomeric interconversion in Cobalt dioxolene 
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We demonstrate that the critical temperature for valence tautomeric interconversion in Cobalt 
dioxolene complexes can be significantly changed when a static electric field is applied to the 
molecule. This is achieved by effectively manipulating the redox potential of the metallic acceptor 
forming the molecule. Importantly our accurate density functional theory calculations demonstrate 
that already a field of 0.1 V/nm, achievable in Stark spectroscopy experiments, can produce a 
change in the critical temperature for the interconversion of 20 K. Our results indicate a new way 
for switching on and off the magnetism in a magnetic molecule. This offers the unique chance of 
controlling magnetism at the atomic scale by electrical means. 



The integration of molecular electronicJi' with 
spintronicfe^ has given birth to the new and fascinat- 
ing field of molecular spintronicn^'^ . This at present en- 
compasses a variety of phenomena, ranging from spin- 
injection and transport in organic semiconductor^ to 
the engineering of magnetic surfaces by means of or- 
ganic molecule^^. A fascinating branch of molecular spin- 
tronics concerns electron transport in single magnetic 
molecule^. Here the aim is that of investigating the 
mutual influence between the spin state of the molecule 
and the electrical current that flows across it. 

The detailed knowledge of the I-V curve of a two- 
terminal devic e m ay enable one to read the spin state 
of the molecul^^H, Much more challenging is the task 
of manipulating such a spin state. Magnetic fields do 
not have enough spatial resolution to address selectively 
single molecules, hence the possibility to use either elec- 
tric fields or currents appears particularly attractive. Re- 
cently it was proposed that spin crossover can be achieved 
by electrical means, simply by exploiting either the firsiP 
or the second order^- Stark effect. The idea is that the 
high-spin (HS) and the low-spin (LS) state of a molecule 
in general have different polarizabilities and different per- 
manent electrical dipoles. Thus the energy Stark shift 
depends on the magnetic state and one can speculate 
that there exists a particular condition where a HS to 
LS crossover is possible. This would have spectacular 
consequences on the molecule I-V characteristic!^. 

Intriguingly there is an entire class of magnetic 
molecules, mainly incorporating Fe(II), known to exhibit 
temperature-induced LS to HS spin-crossover (SCj^. 
Such entropy driven transition is regulated by the rel- 
ative Gibbs free energy between the HS and LS states, 

AG - Ghs - Gls = Ai/ - T^S , (1) 

where AiJ = H^s — Hi^s and AS* = S'hs^'S'ls are respec- 
tively the enthalpy and the entropy variation (AG > 
means that the thermodynamically stable configuration 
is LS). In general in SC molecules AH > but S'hs > 
Shs, so that for temperatures large enough the entropic 
term dominates over the enthalpic one and the molecules 
transit from a LS to a HS configuration. The microscopic 
mechanism for the SC consists in moving an electron from 



a non-bonding to an antibonding state. The transfer pro- 
duces the breathing of the metal ion coordination sphere 
with consequent phonon modes softening. Thus, both the 
spin and the vibrational entropy of the HS state dominate 
over those of the LS ones. Interestingly, it was recently 
demonstrated that the SC is strongly affected by electro- 
static perturbations in the crystal environmenlP^. 




FIG. 1: The Co-dioxolene complex investigated in this work. 
The entire cationic unit of the complex Me2tpa=methil 
derivatives of tris(2-pyridylmethyl)amine, DBCat = 3, 5-di- 
tert-butylcatecolato is presented and the o-dioxolene is en- 
closed in the dashed box. Color code: C=yellow, 0=red, 
Co=grey (large sphere), N=grey (small sphere), H=blue. 

An effect closely resembling the SC is observed in 
Co-dioxolene complexes, wh ich e xhibit an interconver- 
sion between redox isomerJSEH This is called va- 
lence tautomeric interconversion (VTI). The complex 
[Co(Me2tpa) (DBCat)] (PFe) (Me2tpa=methil derivatives 
of tris ( 2-pyrid ylmet h yl) amine , DBCat = 3, 5-di-tert- 
butylcatecolatojfi^E^, is an example of these compounds 
(see Fig. [T]). There are two key ingredients in such a 
class of molecules: the transition metal ion and the o- 
dioxolene group. In general, o-dioxolenes bind to metal 
atoms in three different oxidation states: di-negative cat- 
echolate (Cat), mono-negative semiquinonate (SQ) and 
neutral quinone (Q). The low covalent character of the 
bond between o-dioxolene and Co allows charge and spin 



2 



to be localized so that the molecule maintains well de- 
fined oxidation states and magnetic moments. 

[Co(Me2tpa)(DBCat)](PF6) is found in crystals as ei- 
ther diamagnetic Co(III)-catecholate (LS Co"^+-Cat) at 
low temperature, or as paramegnetic high-spin Co(II)- 
semiquinonate (HS Co^+-SQ) at high-temperature. The 
transition from LS Co'^"'"-Cat to HS Co^+-SQ is driven by 
an intramolecular electron transfer from the catecholate 
to the Co ion. Similarly to standard SC also the VTI is 
entropy driven and can be understood with equation ([ij , 
where LS and HS now correspond respectively to LS 
Co3+-Cat and HS Co2+-SQ. 

Much recent effort has been dedicated to increasing the 
cooperative nature of the VTI in crystals for data stor- 
age applications. Co-dioxolene complexes however might 
also find an interesting place in molecular spintronics. 
In fact the strong interplay between charge transfer and 
spin transition suggests that electric fields might be able 
to profoundly affect the molecule magnetic properties. 
This is demonstrated here. Our density functional the- 
ory (DFT) calculations show that charge transfer with 
consequent spin crossover can be induced by a static elec- 
tric field. Although we find that the critical fields for the 
VTI are rather large (~ 1 V/A), we also show that much 
more modest fields (~ 0.05 V/A) are enough for chang- 
ing the VTI critical temperature by as much as 100 K. 
This intriguingly new effect, which may allow one to fine 
tune electrostatically the magnetic state of a molecule, 
can be readily verified experimentally. 

Achieving an accurate quantitative description of the 
Co-dioxolene electronic structure is a difficult task. One 
in fact has to simultaneously assign the oxidation states 
of both the Co and the dioxolene group, the Co spin 
state and also, in the case of HS Co^^-SQ, the magnetic 
exchange coupling between Co and SQ. Local and semi- 
local approximations to DFT arc inadequate since they 
largely overestimate the electronic coupling between the 
Co 3c? shell and the molecular orbitals of both SQ and 
Cat. This is the same shortfall encountered when describ- 
ing F centers in wide-gap semiconductors and it is rooted 
in the self-interaction erroJ^. In order to overcome this 
limitation we have employed two approaches. The first 
is the atomic self-interaction correction (ASIC) schem^i^l^ 
so far used for electron transport^Q^ and for a broad range 
of problems in solid^^U xhe second is the standard 
B3LYP hybrid exchange correlation functionaPSl. 

Calculations are performed with a development version 
of the SIESTA cod^ implementing the ASIC scheme. 
Norm-conserving Troullier-Martin pseudopotentials are 
employed together with a basis set of double-zeta plus po- 
larization quality. In addiction the GAMES^^ and the 
NWCHElVp^ quantum chemistry packages are used for 
the B3LYP calculations. In this case the basis is 6-31G*. 
We consider the experimental molecular geometries de- 
termined by standard X-ray data analysis for single crys- 
tals (Table ]) without further geometrical optimization 
except for the H atoms. This is the best computational 
strategy since it is at present unclear whether DFT has 



enough accuracy to correctly describe the metal-to-ligand 
bond lengths across the VTI. Note that the same proce- 
dure is often adopted when studying SC compoundd^. 



Isomer 



dco-o (A) 



dco-N (A) 



E (cV) 



LS Co^+-Cat 1.883, 1.898 2.0, 1.947, 2.011, 2.022 
HS Co^+-SQ 2.004, 2.083 2.169, 2.160, 2.1474, 2.151 2.229 
LS Co^+-SQ 1.883, 1.898 2.0, 1.947, 2.011, 2.022 2.789 

TABLE I: Co-O (dco-o) and Co-N (dco-N) bond lengths for 
LS Co^+-Cat and HS Co^+-SQ as determined by X-Ray and 
used in our calculations. The last column reports the B3LYP 
total energy of the given configuration, calculated with re- 
spect to the total energy of LS Co'^^-Cat. LS Co^'''-SQ is a 
fixed spin configuration calculated at the molecular geometry 
of HS Co2+-SQ. 

Table I displays the energy difference between LS Co'^^"- 
Cat, LS Co2+-SQ and the HS Co2+-SQ showing that 
B3LYP qualitatively reproduces the relative order of the 
states. The absolute values of the relative energies appear 
however drastically overestimated with respect to experi- 
mentally available enthalpies (a few hundreds meV)P32ll, 
Such a quantitative disagreement between DFT and ex- 
periments, pathological to many SC compound4^^^'^, 
cannot be attributed to having neglected the vibrational 
zero-point energies, which are much smaller than the elec- 
tronic ones. At present the reasons for these quantitative 
DFT failures are unclear but importantly they affect only 
marginally our main results. 

The ASIC calculated density of states (DOS) of the 
molecule is displayed in Fig. [2] Although in general the 
details of the Kohn-Sham spectrum depend on the spe- 
cific exchange and correlation functional used, in this case 
the relative position of the Go-id shell with the molecu- 
lar orbitals of Cat (SQ) is returned consistently by both 
ASIC and B3LYP. Since, to our knowledge, there are no 
spectroscopic data for Co-dioxolene complexes to com- 
pare with our calculations. Fig. [2] should be interpreted 
as a schematic energy level diagram of the molecule. We 
notice that the self-interaction correction returns the ex- 
pected weak hybridization between the Cat (SQ)-7r* and 
the Co 3d states. The Co "id shell can be clearly rec- 
ognized in the DOS and it is filled according to the for- 
mal oxidation. The symmetry of the highest occupied 
molecular orbital (HOMO) and the lowest unoccupied 
molecular orbital (LUMO) can be appreciated by look- 
ing at charge density isovalue contour plot. For instance 
(Fig. [2| for LS Co3+-Cat the HOMO is the Cat-7r* state 
while the LUMO corresponds to the Co eg states. This 
confirms the good level of the ASIC description. 

After having demonstrated that DFT provides a sat- 
isfactory description of the molecule electronic proper- 
ties we now proceed to investigate the infiuence that an 
electric field, £, has on the VTI. In particular we con- 
sider the situation where £ is applied along the direc- 
tion joining Co with Cat (SQ). In order to reduce the 
computational overheads we replace the N ligands with 
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FIG. 2: Density of states projected over the 3d states of the 
Co and the 2p states of the two O atoms, as calculated with 
ASIC, for a) LS Co^+-Cat b) HS Co2+-SQ c) LS Co'^+SQ 
and d) HS Co'^"*'-SQ with antiferromagnetic coupling between 
the magnetic moment of the Co and the one of the SQ. In the 
top panels we show the isovalue contour plot for the charge 
density (local DOS) for the HOMO (left) and LUMO (right) 
of the LS Co^+-Cat. 



ammonia molecules. This substitution does not bring 
any relevant modification to the relative position of the 
Co-3d and the Cat (SQ)-7r* orbitals. However it has 
an effect on the energetics of the problem, since it de- 
creases the energy difference between HS Co^+-SQ and 
LS Co3+-Cat to 1.591 eV. Figure [s] shows the difference, 
AE = i^HSCo2+SQ - -E'LSCoa+Cat, between the ground 
state energies of LS Co3+-Cat and HS Co2+-SQ {AE > 
means -Bhsco^+sq > -E'LSCoa+Cat)- Such a quantity is 
plotted as a function of the electric field and with respect 
to the £ = situation. The most important observation 
is that AE depends linearly on £. 

We recall here that the energy change, AE, of a quan- 
tum mechanical system under the influence of an external 
electric field is governed by the Stark effect. This reads 




£ (WA) 

FIG. 3: Energy difference between HS Co'^+SQ and LS Co^+- 
Cat, AE, as a function of the electric field £ (the energies 
are calculated with the B3LYP functional). The black line 
(square symbols) corresponds to electric field pointing toward 
the quinone, while the red one (diamond symbols) to the elec- 
tric field pointing toward the Co. The figure also report _Rt 
for the same molecule on right-hand side y-axis. 



AE{£) (X p- £ + ^ '^ij ^i^ij^j where p is the permanent 
electrical dipole and a.y the polarizability tensor. Thus 
the linear dependance found in Fig. [3] indicates that the 
first order Stark effect dominates the molecule response. 
Our DFT calculations confirm such an hypothesis and 
return a finite dipole moment aX £ — Q {\p\ is 16.9 De- 
bye for LS Co3+-Cat and 26.3 Debye for HS Co^+SQ). 
Interestingly the energy decreases when the electric field 
points from the Co to the dioxolene and increases when 
it is along the opposite direction. This means that the 
AE is reduced for the field direction facilitating the VTI 
electron transfer. A linear fit of AE{£) gives us a slope of 
1.12eV/(V/A). We can now extrapolate a critical electric 
field £c = 1.4 V/A at which the states LS Co3+-Cat and 
HS Co^"'"-SQ become energetically degenerate. When £ 
exceeds £c a VTI occurs leading to SC and HS Co^+-SQ 
becomes the new ground state of the molecule. 

The calculated £c unfortunately is quite large, so that 
in reality one can not expect to observe an electric field 
induced VTI. However, it is important to note that such 
a large predicted £c is mainly the consequence of the 
large overestimation of AE{£ = 0), i.e. of the large DFT 
error over the energetics of the system in absence of an 
electric field. A much less sensitive quantity to DFT un- 
certainties is the slope of AE{£), which depends only on 
the molecule electrical dipole and polarizability (accu- 
rately described in DFT). This will allow us to estimate 
with precision the relative change in the VTI crossover 
temperature, Tg, as a function of £ . 

Tq is defined by the condition AG = 0, i.e. it corre- 
sponds to the situation where the Gibbs free energies of 
the two spin configurations become identical. This yields 
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to the equation 



AH 

~AS 



AE 

'as ' 



(2) 



where we have set the pressure to zero since we are con- 
sidering the molecule "in vacuum" , and we have approx- 
imated the enthalpy at Tc with the zero-temperature to- 
tal electronic energy. In writing the equation ^ we have 
also assumed that the entropy variation, AS, does not 
depend on the electric field, i.e AS{£) « AS{0). This 
is consistent with our computational strategy of not re- 
laxing the atomic coordinates. We have verified in a 
few cases that the atomic relaxation under bias is much 
smaller then that due to the VTI, so that our assumption 
is well grounded in our own calculations. 

Then, some algebra shows that the ratio between the 
critical temperature calculated in an electric field, Tc{£), 
with that in zero field, Tc(0), simply writes 



i?T — 



Tc{£) ^ AE{£) 
Tc(0) AE{Q) 



(3) 



Importantly is expressed in terms of DFT total en- 
ergies only, and in Fig. [3] we show that changes of Tq 
of the order of 5% can be obtained already at the much 
more modest field of 0.07 eV/A. 

The uncertainty over such a value nevertheless remains 
quite large since the calculation still involves the pre- 



cise determination of AE{Q). A better estimate can be 
obtained by eliminating AE{0) completely and just use 
information about the slope of AE{£). This can be ob- 
tained by re-writing the equation ([2| as 

[Tc(0) - Tc(f )] A5(0) = AE{Q) - AE{£) . (4) 

We now take the experimental estimat^^H for the en- 
tropy variation AiS'(O) « 0.52 nieV/K, and conclude that 
achievable electric fields can indeed induce large varia- 
tions in Tq. In fact an electric field of 0.1 V/nm, typical 
of Stark spectroscopjPsl^ produces already a change in 
Tc of about 21 K, and a field of 0.5 V/nm, obtainable in 
scanning tunnel microscopy experiments can drift Tq by 
as much as 100 K. 

In conclusion we have demonstrated that the critical 
temperature for the VTI in Co-dioxolene complexes can 
be drastically modified by a static electric field. Intrigu- 
ingly, this represents a physical way to effectively modu- 
late the redox potential of the metal acceptor, alternative 
to the chemical strategy of changing the radical groups 
of the molecule--. Such an external control of the mag- 
netism of valence tautomeric compounds may open new 
avenues to the growing field of molecular spintronics. 
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